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Abstract — In diffusion-based molecular communications, mes- 
sages can be conveyed via the variation in the concentration of 
molecules in the medium. In this paper, we intend to analyze 
the achievable capacity in transmission of information from one 
node to another in a diffusion channel. We observe that because 
of the molecular diffusion in the medium, the channel possesses 
memory. We then model the memory of the channel by a two-step 
Markov chain and obtain the equations describing the capacity 
of the diffusion channel. By performing a numerical analysis, we 
obtain the maximum achievable rate for different levels of the 
transmitter power, i.e., the molecule production rate. 

I. Introduction 

Our interest to study capacity of molecular diffusion chan- 
nels is twofold. First, there has been numerous evidence of the 
existence of different forms of communication in nature. Com- 
munication enables single cells to process sensory information 
about their environment (in a way similar to neural networks) 
and evaluate and react to chemical stimuli. At the microor- 
ganism scale, molecular signals are used for communication 
and control among cells. For example, one of the well-known 
communication primitives among cells is the phenomenon 
called Quorum Sensing. Quorum sensing is a cell-to-cell 
communication process in which bacteria use the production 
and detection of molecules to monitor the population density 
of the bacteria. Quorum sensing allows bacteria to synchronize 
the behavior of the group, and thus act as a unit fTl-f4^. 
The communication can also appear in other forms such as 
calcium signaling. Cells absorb calcium molecules in response 
to various stimuli that open/close particular channels on the 
cell membrane. The molecular information in the variation of 
the calcium ions concentration is propagated both inside and 
outside the cell, causing a variation in the electrical charge of 
the cell membrane and, subsequently, the transduction of the 
information into an electrical signal. Due to various limitations 
such as size (size of a single cell or microorganism) and energy 
in small scales, the dominant form of communication in such 
scales is via molecular signals, which is fundamentally differ- 
ent from conventional electromagnetic-based communication. 

Secondly, recent advances in bio-nano technology has mo- 
tivated research on designing nanoscale devices to perform 
tasks similar to their biological counterparts [5 |. There is a 
large number of applications that nano-devices could apply 
to. One may envision molecular based networks built using 
these nanoscale devices that can be deployed over or inside 
the human body to monitor glucose, sodium, and cholesterol 
levels, to detect the presence of different infectious agents, or 
to identify specific types of cancer Such networks will also 



enable new smart drug administrative systems to release spe- 
cific drugs inside the body with great accuracy and in a timely 
manner To enable all such applications, the communication 
among nano-devices is the key. 

The above promising outlook has inspired development 
of new theoretical frameworks for molecular communication 
such as 161, Q. These studies by Eckford, et. al. mostly 
focused on molecular communication, with information con- 
veyed by the time of the release of molecules while the prop- 
agation between the transmitter and the receiver is governed 
by Brownian motion (without drift fS] and with drift [7], 
||9l). An alternative molecular communication system is the 
one governed by diffusion process in the medium, called 
diffusion channel. Arguably, the most dominating form of 
the communication at the micrometer scale is diffusion based 
molecular communication, i.e., embedding the information in 
the alteration of the concentration of the molecules and rely 
on diffusion to transfer the information to the destination. The 
communication among bacteria, calcium signaling, pathogen 
localization functionality of the immune system and many 
others can be reduced to diffusion based molecular communi- 
cation. Another body of work in this field involves the study 
of the Quorum Sensing as a network and mapping the Quorum 
Sensing to consensus problem under diffusion-based molecular 
communication ifTOl . 

In this work, we study the diffusion channel. We consider 
a communication scenario consisted of a transmitter and a 
receiver communicating via the diffusion channel, as shown 
in Fig. [U In this context, we can think of the transmitter as 
a nano-device or a bio-engineered cell capable of emitting 
molecules and releasing them into the medium to change 
their concentration according to information bits. Similarly, 
the receiver is capable of absorbing the molecules or chemical 
signals, for example, by using ligand-receptor binding which 
is a transmembrane receptor protein on a receiving cell [TT] , 
||T2|. This ligand-receptor interaction creates peculiarities such 
as non-Gaussian noise. The receiver has a large number of 
binding places using which it can estimate the concentration 
by averaging over all binding places. The diffusion process 
has a profound impact on transmission of information that 
makes the diffusion channel very different from the classical 
models developed for electromagnetic-based communication. 
The molecules that are produced by the transmitter stay 
in the medium and affect the later transmissions. Our goal 
is to compute the maximum achievable rate of information 
exchange (hereafter referred to as capacity) in such a diffusion 
based molecular communication channel. Analogous to the 
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Fig. 1. Scheme of the model for molecular communication 

early studies of communication systems, we avoid some of 
the peculiarities of the system in our first step and we focus 
on studying the simpler and more practical discrete noiseless 
systems. According to Shannon, in a discrete communication 
system, a sequence of choices from a finite set of elementary 
symbols S can be transmitted from transmitter T to receiver R. 
Each symbol Si is assumed to have duration Ts- . As it it will 
be discussed later, it is important to notice that the symbols 
have different durations imposed by the diffusion process. 

Therefore, to compute the capacity, we first characterize, 
the duration of each symbol by a careful study of diffusion. 

In Seciniwe review some results about diffusion process. In 
SecHUwe introduce the communication model and investigate 
the problem under study in more detail. Our main results are 
discussed in Sec. IIVI 

II. Background 

The first step is to characterize the temporal and spatial 
variations of molecules in the channel which follows the 
general diffusion equations. According to Pick's second law 
of diffusion the concentration of molecules c{x, t) at position 
X at time t is computed using the molecule production rate 
r{x,t), as follows: 



dc(x, t) 
dt 



= D\7^c{x,t) +r{x,t), 



(1) 



Here, x is the distance of any point in the environment from the 
source and D is the diffusion coefficient of the medium. The 
impulse response of ([Tjl, corresponding to r{x,t) — S{x)6{t), 
is the Green's function gd{x, t) whose expression is as follows: 
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This impulse response is given for the 2-D medium but 
can be extended to a general 3 -dimensional case, using the 
observation that n-dimensional diffusion is equivalent to n 
separate (simultaneous) 1-D diffusions. Since the diffusion 
equation is a linear equation, the solution to ([Tjl for an 
arbitrary input r(x,i), denoted by c*{x,t), can be obtained 
using 

c*{x,t) = gd{x,t) ®r(x,t), (3) 



where ® denotes a multi-dimensional convolution operation 
on X and t. 

In our setup, we assume that there is only one transmitter 
emitting molecules. Therefore, in 0, we have r{x,t) — 
F{t)S{x), where F{t) is the input signal. Hence, we will have 
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(4) 

This response is valid for open free media in which the 
only boundary conditions are at the transmitter Note that in 
this model, we do not consider the delay due to the travel 
time of molecules between the transmitter and the receiver 
That is we assume molecules reach the receivers instantly. 
This assumption, however, does not affect our analysis of the 
capacity because it only shifts the time that molecules are 
arrived at the receiver 

III. Model 

We consider the case of a single transmitter and a single 
receiver in a 2-D environment. This model can be extended 
to a network of transmitters and receivers in which each 
transmitter intends to communicate with a specific receiver 
when interference among these communications is negligible. 
We consider a discrete-time model for the communication in 
which an arbitrary string of binary information is given to the 
transmitter This string is to be conveyed to the receiver at 
distance r from the transmitter 

The transmitter encodes the information bits into the proper 
concentration of molecules and releases them for some specific 
durations, discussed later These molecules diffuse in the 
medium. The receiver compares the concentration of received 
molecules with a specific threshold and decides whether the 
concentration in the environment is low or high. This threshold 
depends on the characteristics of the environment. It can nei- 
ther be arbitrarily small because of the noise and interference 
in the environment nor be arbitrarily large because of the 
limited capability of the transmitter for molecule production. 
In the latter case, either the concentration may not reach the 
threshold or it would take a long time. 

The main difference here with a typical paradigm of com- 
munication is the memory which essentially exists in the diffu- 
sion channel and influences the communication dramatically. 
The concentration of molecules cannot change momentarily 
and unlike the typical binary symmetric channel, previously 
sent symbols affect the new ones. For example, suppose the 
concentration of molecules in the channel is in the high state 
and the transmitter intends to send "0" through the channel. 
However, molecules from the previous state still exist in the 
environment and cannot be removed suddenly. Hence, the 
transmitter has to wait a specific amount of time to send the 
new bit. Note that we do not assume existence of negative rates 
because it is unrealistic to assume as such for the molecular 
communication. We also do not consider any noise to simplify 
the study. 

Based on the previous discussion, we observe that each bit 
depends only on the last bit that was transmitted and it is 




determinant equation \M\ = 0, i.e.. 



Fig. 2. Two states of the channel and possible inputs in each state. 

independent of the other previously transmitted bits. Therefore, 
a Markov chain can be used to model the communication in 
the diffusion channel. In such a model, the state of the channel 
is set by the last transmitted bit: the channel state is "H" if the 
last transmitted bit was "1", and it is "L" otherwise. Hence, 
there are two states in the channel, as shown in Fig. |2] In the 
low state, we can either stay in this state by sending "0" or 
go to high state by sending "1". In the high state, however, 
for sending the bit "1" again, it is not efficient to wait for the 
channel to reset to "L" state. Instead, the transmitter can send 
"1" which is shorter in time duration and needs less molecules 
or equivalently less energy. To send "0" when channel is 
at state "H", it suffices that the transmitter waits a specific 
amount of time such that the channel clears itself from the 
molecules; returning to "L" state. 

The above discussion implies that depending on the channel 
state, duration of bit "0" is different. Likewise, the duration for 
which the transmitter emits molecules for bit "1" depends on 
the channel state. Therefore, effectively, we can think of four 
different symbols sent by the transmitter, namely, so ="00", 
si ="01"', S2 ="10"', and S3 ="11"'. The first bit in each 
symbol indicates the previous bit which was already sent (and 
hence the channel state), and the second bit is the one to 
be transmitted next. Assuming a binary symmetric source, 
all these four symbols are equiprobable. However, a different 
transmission time is associated with each, as we explained. 

Since our model involves symbols with various time dura- 
tion, we adapt the method by Shannon in ifTSl to compute the 
capacity of a discrete noiseless channel as 



C= lim i^^. 



(5) 



where N{T) is the number of allowed blocks of duration T. 

(s) 

Let ' be the duration of s-th symbol which is allowable in 
state i and leads to state j. Then Nj (T), the number of blocks 
of length T ending in state j is given by 



iv,(T) = ^iv,;(r-4^') 



(6) 



The asymptotic solution for these difference equations 
would be in the form Nj = AjW'^ where Aj is a con- 
stant 1 13 1. By substituting this form in the possible number 
of blocks would be equal to '^jAjW'^ where W can be 
found by solving the equation ^ below. We define a matrix 
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Here, 5ij = 1 if i = j and is zero otherwise. Hence, 
following (|5]l, the channel capacity C is obtained as 
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Rewriting O for the diffusion channel described above results 
in 



= 0. 



(8) 



Here, where s e {00,01,10,11} is the duration of the 
symbols. This leads to 

By solving Q with respect to Ts for s G 00,01, 10, 11, we 
obtain the capacity of the channel in bits per second. In the 
next section, we use the diffusion channel property to solve 
for T,. 

IV. Main Results 

Consider the scenario in which for symbol s, the transmitter 
sends a pulse-shaped rate of molecules with amplitude Fg to 
the receiver for duration T,. It is clear that Fg is zero for 
symbols si = 00 and S2 = 10 because the transmitter does 
not need to emit any molecules. Instead, it must wait specific 
amounts of time. We assume Fqi (for si) to be equal to F, the 
maximum rate that the transmitter can produce. For S3 = 11, 
there is no need for the transmitter to send as many molecules 
as for si. Hence, we assume Fn = aF. Thus a fraction a of 
the maximum rate is allocated to S3. We denote by Cg{r,t) 
the response at the receiver, i.e. x = r, at time t due to 
transmission of symbol s. Using (|4|, we have 
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5ij]\ W is the largest real root of the 



(10) 

The diffusion equations of the channel do not set any limits 
on Too. It depends on the receiver sensitivity, the distance 
r between the transmitter and receiver and the diffusion 
coefficient D. Specifically, Tqo can be considered as the time 
sensitivity that the receiver can sense the medium. Hence, Too 
is the fundamental time. The other T, for s G {01, 10, 11} can 
be considered as a multiple of Too. Based on the Markov chain 
model, we derive the equations for Tg for s G {01,10,11} 
corresponding to different symbols. 

Assume channel is in the low state. The transmitter may 
send bit "0", i.e sq = 00, or change the channel state to high, 
i.e. si = 01. In the former case, the transmitter must wait 
for an interval of Too but in the latter, it must emit molecules 



into the channel. Hence, for the si = 01, we will have (for 

0<t< Toi) 

c*{r,t)^coi{r,t). (11) 

On the other hand, when the channel is in the high state, 
we have to take into the account the effect of the existing 
molecules in the channel, i.e. the memory of the channel. We 
can replace the molecules in the channel by the source that 
produces this concentration before the transmission of the next 
symbol and is zero afterwards. Based on our model, we can 
capture the memory of the channel by incorporating the effect 
of the symbol si on the next symbol. In the case that the 
transmitter wants to send "0" afterwards, i.e. S2 = 10, it must 
wait for the channel to become low. Hence, we have (for < 
t < Tio) 

c*{r,t) = CQi{r,t + Toi) + cio(r,t) = CQi{r,t + Tqi). (12) 

where the first term corresponds to the effect of the memory. 
The second term cio{r,t) is equal to zero which implies that 
the transmitter does not need to send any molecules but has 
to wait an appropriate amount of time. Finally in the case that 
the transmitter wants to send "1" again, i.e. S3 — 11, we have 
(for < t < Tu) 



'{r,t) ^ Coi{r,t + T01) + cn{r,t). 



(13) 



Here the first term in (ITJt corresponds to the effect of the 
channel memory and the second term is due to the emission 
of molecules for this symbol. 

In the absence of noise, any intended concentration can 
reach the receiver without any interference from the medium. 
Therefore, the only limiting factor for capacity is the memory 
of the channel imposed by diffusion transmission. Let S denote 
the concentration sensitivity of the receiver This implies 
that the receiver cannot differentiate between the levels of 
concentration that differ by less than S. We map the low 
concentration on S and high concentration on 25. We want 
the concentration at the receiver to be equal to one of these 
values at the end of each interval, i.e. t ~ Tg. Using (fTTl l. (fT2l) . 
and ( fT3]) . we have 

coi{r,Toi) = 2S, 

coi{r,Toi+Tw) = S, . (14) 

coi('','7oi + Tu) + cii(r,Tii) = 25" 

Note that in these equations the time delay between the 
transmitter and the receiver is not considered. We can write the 
solution for (fT4] i based on the function Ei(a;) defined as I14ll 
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Then expressions in ( fT4l i can be written as 
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Fig. 4. Impulse Response of the Diffusion Channel. 

Since finding closed form formula for (fT&i is not possible, 
in the following we solve for ( fTSI l numerically to obtain Tqi, 
Tio, and Tu in terms of F. The parameter F indicates the 
maximum power used by the transmitter which controls the 
capacity. We solve ( fTSI ) based on the normalized parameter 

~ 4ttDS ^'^'^ obtain Tg for s £ {01, 10, 11} as multiples of 
jj^. As discussed before, Tqo depends on r and D. Therefore, 
we may assume Too = where fc is a constant. We 

arbitrarily choose fc = 1, although our analysis holds for any 
other values of fc. 

It can be shown that Tu is smaller, by orders of magnitude, 
than both Toi and Tio. Intuitively, this is due to the fact 
that producing a high state needs much more molecules than 
maintaining it. Further, as indicated in the following, Toi and 
Tio are in the same order of magnitude for typical values of 
F. Hence, for simplicity, we assume Tu to be equal to Tqo, 
and perform our numerical analysis only on Toi and Tio. The 
numerical results for Tqi — and Tio — ^ are shown in 



Too 



Fig. |3(a)| and |3(b)| respectively. 

Fig. [3] shows the plots of transition time versus the normal- 
ized maximum rate of the molecule production F. We observe 
that the transition time from low to high Toi is monotonically 
decreasing with respect to F, and approaches to infinity when 
F is close to zero. In Fig . |3(b)| the normalized transition time 
from high to low state Tio is presented. It has a minimum 
around F = 2. This behavior can be explained by the diffusion 
equation as following. Using (|2]i, the impulse response of 
the diffusion channel is shown in Fig. |4] As we see in this 
plot, the impulse response increases to a maximum and then 
falls lower. To obtain the response to the transmitter, a pulsed 
shape signal in time has to be convolved with the impulse 
response. When we increase the F (or equivalently decrease 
Toi), at some point the threshold 25 is reached before the 
maximum of impulse response. Since the tail of the impulse 
response decreases slowly, in order to send the symbol "10", 
the transmitter has to wait until the pulse is slided to the right 
side of the maximum. When F is smaller and Tqi is larger, 
we will not face this problem. Therefore, using Fig. |3(b)| we 
conclude that increasing the rate of production F does not 
always result in smaller Tio. Although larger F results in 
smaller Toi, it makes Tio to grow larger, and hence has a 
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Fig. 5. Capacity vs Production Rate. 



negative effect on the capacity. Hence, the optimum value for 
the production rate can be found by trade-off between Toi and 
Tio. 

Using ^ and the assumption that Tn = Tqo, the capacity 
of the channel can then be obtained by solving for W in the 
following equations. 
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and using it in C = logVF. The capacity, logVF, is in "bits 
per an interval of Too". We note that maximum capacity 
is achieved when Tqi + Tio is minimized. The plot of the 
resulting capacity is shown in Fig. |5] with respect to F. As 
shown in the figure, the maximum capacity is approximately 
0.5 (bit/Too) which is achieved when choosing F = 3.9 (units) 
and corresponds to the minimum value of Tqi + Tig which 
is approximately 7 (time units). The left tail of the graph 
converges to zero which corresponds to the no production 
rate. The right tail of the graph converges to zero as well 
but at a slower pace. In this case, symbols "00", "11", and 
"01" can be transmitted assuming the channel is in the proper 
state. However, since Tip becomes very large, the transition 
between the channel states takes a long time, leading to zero 
capacity. ^ Conclusion 

In this paper, we used the discrete noiseless channel model 
to obtain the capacity of the diffusion channel. The inherent 



memory in the diffusion channel distinguishes it from typical 
wireless channels, making the capacity analysis more challeng- 
ing. For a binary input, we considered four different symbols 
to take into account the effect of the channel state on the 
next transmitted bit. With a numerical analysis, we showed 
that beyond a threshold, increasing the level of power has a 
negative effect on the capacity, which was explained by the 
physical characteristics of the diffusion channel. 
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